NORDITA - 94/11 N,P 



TWO- AND THREE-POINT FUNCTIONS 
IN THE EXTENDED NJL MODEL 



Johan Bijnens" and Joaquim Prades"' 6 

a NORDITA, Blegdamsvej 17, 
DK-2100 Copenhagen 0, Denmark 

6 Niels Bohr Institute, Blegdamsvej 17, 
DK-2100 Copenhagen 0, Denmark 



Abstract 

The two-point functions in generalized Nambu-Jona-Lasinio models are 
calculated to all orders in momenta and quark masses to leading order in 
1/N C . The use of Ward identities and the heat-kernel expansion allows 
for a large degree of regularization independence. We also show how this 
approach works to the same order for three-point functions on the ex- 
ample of the vector-pseudoscalar-pseudoscalar three-point function. The 
inclusion of the chiral anomaly effects at this level is shown by calculat- 
ing the pseudoscalar-vector-vector three-point function to the same order. 
Finally we comment on how (vector-)meson-dominance comes out in the 
presence of explicit chiral symmetry breaking in both the anomalous and 
the non-anomalous sectors. 
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1 Introduction 



The Extended Nambu-Jona-Lasinio model (EN JL) [0], ^ has already a long his- 
tory. For some recent reviews see [[| and references therein. Generally a good 
agreement with low-energy hadronic phenomenology has been found. However its 
main drawback is the lack of confinement. In ref. Q a large number of relations 
between the observables was found which were valid in a large class of ENJL-like 
models. In ref. || this type of relations was generalized to two-point functions 
and to all orders in the momenta in the chiral limit. Various numerical results 
obtained in ref. || thus obtained a larger range of validity. In this work we shall 
extend this type of analysis to three-point functions and two-point functions be- 
yond the chiral limit. The first one, the vector-pseudoscalar-pseudoscalar can 
already be found in |?j and we have included it to show explicitly the use of one- 
loop Ward identities to simplify the calculation. A similar approach can be found 
in ref.fl which we received when the analytic part of this work was essentially 
finished. They have a less general treatment of regularization dependence than 
is done here and only treat the case with equal current-quark masses. With the 
same definitions our results for the two-point functions agree with theirs. Our 
main aim is to apply this procedure to the case of the pseudoscalar- vector- vector 
three-point function. Here we both illustrate our prescription for the consistent 
treatment of the chiral anomaly in this model by imposing the QCD anomalous 
Ward identities M. The latter do imply the use of consistent one-loop ENJL 
anomalous Ward identities. We find that the duality between the Vector-Meson- 
Dominance (VMD) picture for the slope of the anomalous 7r°77 form factor and 
the quark-loop one is much worse that the one found for the pion electromagnetic 
form factor and a more refined model (like ENJL cut-off like models) is necessary 
to reconcile both approaches. 

At this point we would like to add some comments about the anomaly in the 



ENJL model. In [ ID | it has been argued that the anomaly can not be consistently 
reproduced in this type of models. While we agree that in general in these models 
no simple definition of the anomaly is possible, we believe, as discussed in ||, that 
if one wants to use these models as a low-energy approximation to QCD there is 
a unique prescription of how to do this. Other uncertainties in the regularization 
scheme of the anomalous sector will be suppressed by powers of the cut-off A x 
used in the ENJL model. 

The paper is organized as follows. In section |2] we give a short overview of 
the ENJL model. Here we also discuss the dependence of the constituent quark 
mass on the current quark mass and make some remarks about the definition of 
quarks versus the QCD ones. In section [3] we extend the analysis of ref. |j| to 
the case with nonzero current quark masses and both masses that play a role in 
the two-point function are allowed to be different. We compare the results with 
Chiral Perturbation Theory (xPT) and show numerical results for some of the 
two-point functions. We also discuss the method used shortly and give the new 
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identities that the two-point functions at one-loop and the full resummed ones 
have to satisfy. Their derivation is rather technical and has been given explicitly 
in appendix [B]. The main difference with ref. is that now there is also non- 
trivial mixing in the scalar sector. In subsection Q.8 we discuss in detail the 
Weinberg Sum Rules (WSR). It is found that here the high energy behaviour of 
this class of ENJL-like models is too strongly suppressed. 

Then we come to the derivation of the three-point functions here in the next 
section [|. In subsection |0| we discuss the vector-pseudoscalar-pseudoscalar 
three-point function paying attention to the Ward identities in its calculation. 
In subsection |12] we do the same for the pseudoscalar- vector- vector three-point 
function. Here we explain how one needs to treat the anomalous part of the Ward 
identities to get a consistent result. Section |5| treats the appearance of Vector 
Meson Dominance like features of two- and three-point functions in this class of 
models. We briefly discuss two-point functions and particularly the transverse 
vector two-point function in the first subsection. Here the origin of the large shift 
in the slope compared to My(0) of ref. M is explained. In subsection E]2 



we 



discuss the VMD behaviour of the first three-point function and give a discussion 



of the KSRF identity [[LI]] in this ENJL model. In the last subsection we treat 
the vector-pseudoscalar-pseudoscalar three-point function similarly. In section ^ 
we summarize our results. 

The appendices contain the definition of our regularization procedure, the 
derivation of the Ward identities and explicit expressions for the one-loop func- 
tions we need. 



2 Short description of the ENJL model and its 
connection with QCD 

The QCD Lagrangian is given by 

£qcd = £qcd — ^G^G^ , 

£qcd = Q {vf ~ iv » ~ ia »l5 ~ -(M + s- ijry 5 )} q . 

(2.1) 

Here summation over colour degrees of freedom is understood and we have used 
the following short-hand notations: q = (u, d, s) ; G^ is the gluon field in the 
fundamental SU(N C ) (N c =number of colours) representation; G^ v is the gluon 
field strength tensor in the adjoint SU(N C ) representation; t> M , a^, s and p are 
external vector, axial-vector, scalar and pseudoscalar field matrix sources; M. is 
the quark-mass matrix. 

All indications are that in the purely gluonic sector there is a mass-gap. There- 
fore there seems to be a kind of cut-off mass in the gluon propagator (see the 
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discussion in ref. |T2|). Alternatively one can think of integrating out the high- 
frequency (higher than A x , a cut-off of the order of the spontaneous symmetry 
breaking scale) gluon and quark degrees of freedom and then expand the result- 
ing effective action in terms of local fields. We then stop this expansion after the 
dimension six terms. This leads to the following effective action at leading order 
in the 1/N C expansion 

£qcd — > £qcd + £njl + ^njl + C (l/ A£) , 

with 41 = Sj ^r A T,(?W L )(^' R ) 



and 



(2.2) 



Where i, j are flavour indices and St^z, = (1/2) (1 ±75) ^. The couplings Gs 
and Gy are dimensionless and 0(1) in the l/iV c expansion and summation over 
colours between brackets is understood. The Lagrangian Cqq B includes only low- 
frequency modes of quark and gluon fields. The remaining gluon fields can be 
assumed to be fully absorbed in the coefficients of the local quark field operators 
or alternatively also described by vacuum expectation values of gluonic operators 
(see the discussions in refs. 0, 3). In the mean-field approximation these £njl, V ' A 
above are equivalent to a constituent chiral quark- mass term ||13|| . 

This model has the same symmetry structure as the QCD action at leading 
order in 1/N C [14] (notice that the U(1)a problem is absent at this order [15]). 



(For explicit symmetry properties under SU(3)l x SU(3)r of the fields in this 
model see reference [|J.) We can self-consistently solve the Schwinger- Dyson 
equation for the fermion propagator in terms of the bare propagator and a one- 
loop diagram. In the case where the current quark masses are set to zero this 
equation allows for two solutions for Gg > 1, one with constituent quark mass 
M = and the other with M 7^ and the model shows spontaneous chiral 
symmetry breaking. In the presence of explicit chiral symmetry breaking only 
the second solution is allowed. We shall allow for nonzero current quark masses, 
M. = diag (m u , m d , m s ) and all different. In the leading l/N c limit the solution 
of the Schwinger-Dyson equation is a flavour diagonal matrix for the constituent 
quark masses with elements M u ^ s - The gap-equation now becomes 

Mi = mi -g s (0\:q iqi : |0) , (2.3) 

(0| : m : |0> = (q^i) = -NAM, f ^ 

N, 



16n 



(2vr)V -M? 
2 4M?r(-l, ei ) , (2.4) 
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m = MeV 
m 10 MeV 




GS 

Figure 1: Plot of the dependence of the constituent quark mass Mj as a function 
of Gs for several values of m; 



» s wi ■ < 2 - 5 > 

Therefore, in this model the scalar quark-antiquark one-point function (quark 
condensate) obtains a non-trivial nonzero value. The dependence on the current 
quark- mass is somewhat obscured in eq. Q2.4Q . We use here a cut-off in proper 
time as the regulator. See appendix |A] for its definition. The quantity appearing 
in ( |2.4| ) is Mf/A^. In figure [I] we have plotted the dependence of Mj on Gs for 
various values of m« and A x = 1.160 GeV. It can be seen that the value of 
Mi for small rrii converges smoothly towards the value in the chiral limit for the 
spontaneously broken phase. This is an indication that an expansion in the quark 
masses as Chiral Perturbation Theory assumes for QCD is also valid in this model. 
However, it can also be seen that the validity of this expansion breaks down 
quickly and for mj ~ 200 MeV we already have 2Mj ~ A x . We note that the ratio 
of vacuum expectation values for light quark flavours increases with increasing 
current quark mass at p 2 = in this model and starts to saturate for mj > 200 
MeV. In standard this ratio is taken to be 1 at lowest order and its behaviour 
with the current quark mass is governed (at 0(p 4 )) by the following combination 



of coupling constants 2L$ + H2 [I6j in the large N c limit. The 0{p ) %PT coupling 



constants |T6j are calculated at leading order in 1/N C and in the chiral limit in 



the ENJL model Wp. The analytical result for this combination of couplings 



lr Tlie analytical expression for if 2 in that reference is correct. The tables contain a numerical 
error. For example the value of H2 for the parameters of fit 1 in ref. B is 1.4 • 10 -3 
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constants there was confirmed in ref. H from a calculation of the scalar two-point 
function in the chiral limit and in the large N c limit to all orders in momenta. 
This result was obtained by requiring the relevant Ward identities to all orders 
in momenta. The value for the combination 2L§ + found there corresponds to 
a big increase of the quark condensate with the current quark mass at this order. 
We want to emphasize here that the exact identification of the quark condensate 
in eq. (|2.4j ) which is regularization dependent and especially its dependence on 
the current quark mass with the one used in %PT or QCD Sum Rules is by no 
means straightforward. The differences between both can be traced back in the 
scalar sector of the model and in particular in the quadratically regularization 



dependent 0{p ) coupling constant H 2 [fU|. This high-energy constant is related 



to the details of the integration of the QCD high-frequency modes to obtain the 
Lagrangian in eq. Q2.2|) . However there are indications that the QCD light quark 
condensates indeed increase with the current quark mass. In general, there is 
some uncertainty in the definition of the quark fields in ENJL models versus the 
QCD ones. This depends on the details on how the ENJL model originates from 
QCD. 



3 Two-point functions in the presence of cur- 
rent quark masses 

This section is a generalization of the results in ref. [||] to the case of nonzero 
current quark masses. These two-point functions were studied before in || but 
there they were discussed as quark form factors. What is new here is that the 
explicit dependence on the regularization scheme has been put into two arbitrary 
functions, namely, Tly^ + Tly^ and II^j (see this section below for definitions). 
This also shows that these results are valid in a class of models where the one- 
loop (see further for the definition of this) result can be expanded in a heat-kernel 
expansion using the same basic quantities E and R^ v as used here. This includes 
the ENJL model with low-energy gluons described by background expectation 
values. We have not included this case in our numerical results for the explicit 
one-loop expressions. For the equal mass case the relevant one loop formulas can 
be found in ref. 0. 



3.1 Definition of the two-point functions 

We shall discuss two-point functions of the vector, axial-vector, scalar and pseu- 
doscalar quark currents with the following definitions, 

V;\x) = ft(x) 7 ^(x), (3.1) 
A^(x) = qi(x) W qj(x) , (3.2) 
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S ij (x) = -qi{x) qj (x), (3.3) 

P^{x) = q^i^x), . (3.4) 

The indices i,j are flavour indices and run over u,d,s. The two-point functions 
themselves are defined as 

KMm = ^/d 4 xe^<0|T(^'(x)^(0))|0>, (3.5) 

n£,(ff)«« = ^/d 4 xe^<0|T(^(x)^(0))|0>, (3.6) 

nJ(g) iiW = i / d 4 xe 1 ^ < 0|T (v™(z)S w (0)) |0 > , (3.7) 

K{q)w = 1 J d ' xet9 X < °I T 1° > > ( 3 - 8 ) 

n 5 (g) ijW = i J d'xe 1 ^ < 0|T (S l > (x) S kl (0)) |0 > , (3.9) 

ll p {q) m = % J d 4 xe iq x < 0|T (p»{x) P kl (0)) |0 > . (3.10) 



In the leading order in the number of colours these are all proportional to Sijki = 
duSjk, with 5u the Kronecker delta. Using Lorentz-invariance these functions can 
then be expressed as follows 

u lv{(l)ijki = {{q^qu-q 2 g^ v )^v\Q 2 )ij + %qu^v\Q 2 )ij} 5 ijkU (3.11) 

n%{q)w = {(q^-q%u)iiT(Q\ + Q^iif(Q 2 h}s ijkl , (3.12) 



nfa)w = q^s(Q 2 h s m, (3.13) 

nj(g)ii« = iq^p{Q%AikU (3.14) 

n 5 ( g ) ijM = n s (g 2 )i/ iiW , (3.15) 

U p (q) ijkl = Iip{Q 2 ) l3 5 l]kl . (3.16) 



Here Q 2 = —q 2 . We shall discuss the Weinberg Sum Rules and numerical results 
for the two-point functions only in the Euclidean domain, i.e. Q 2 positive. Us- 
ing Bose symmetry on the definitions of the two-point functions it follows that 
W(Q%, Ii ( v\Q 2 )i V ni 0) (Q 2 )^, n?(Q 2 ) li; U S (Q% and U M (Q% are all sym- 
metric in the flavour indices i and j. The remaining ones need the Ward-identities 
to prove their flavour structure. From the identities in the appendix [B] it follows 
that Ug I (Q 2 )ij is also symmetric in while IT^Q 2 )^ is anti-symmetric. 

3.2 Lowest order results in Chiral Perturbation Theory 

From Chiral Perturbation Theory to order p 4 in the expansion we obtain the 
following low energy results for the two-point functions. The orders mentioned 
behind are the orders in Chiral Perturbation Theory that are neglected. 

n { v\Q 2 h = -^2H 1 + L 10 ) + O(p 6 ), (3.17) 
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li 0) (Q% = (9(p 6 ), (3.18) 

^a\Q% = 2 ^-i(2H 1 -L 10 )+O(p e ), (3.19) 

nf(Q% = 2f> (3.20) 

n^(Q% = 0(/), (3.2i) 

n£(Q 2 ) y = 4^% + b 6 )> (3-22) 



n s (Q% = 85 2 (2L 8 ±# 2 ) + O(p 6 ), (3.23) 
Kp(Q 2 h = 2 ^!% + 8/j Q 2 (-2L 8 + H 2 ) + 0(p 6 ) . (3.24) 

With rriij the mass of the lightest pseudoscalar meson with flavour structure ij. 
These are obtained in the leading 1/N C approximation so loop-effects are not 
needed. Notice that these expressions are valid to chiral order p A . From a term 
of the form ^{D^xD^X 1 } there are contributions of order (mj — rrij ) 2 /Q 2 to the 
vector two-point function Tly (Q 2 )ij and of order (m^ — rrij) to the mixed scalar 
vector function Ilcf (Q 2 )r/. 

The functions n^, lip and lip get their leading behaviour from the pseu- 
doscalar Goldstone pole. In addition IT^ and IT^ contain a kinematical pole at 
Q 2 = 0. The residue of the physical pole is proportional to the decay constant 
fij for the relevant meson, (for the ud ones, f u d ~ f n ~ 92.5 MeV). In x?T, the 
constant Bq is related to the vacuum expectation value in the chiral limit. In 
the large N c limit and away from the chiral limit there are corrections due to the 
terms proportional to combination of 0(p A ) couplings 2L 8 + H 2 [IB|. 



< 0| : W : |0 > l9=uAs = -f%B (l + 0(p 4 )) . (3.25) 

The vacuum expectation value here, < 0| : : |0 >, is the one used in x?T in 
the chiral limit and fo is the pseudoscalar meson decay constant in the chiral limit. 
The constants L 8 , L 10 , Hi and H 2 are coupling constants of the 0(p 4 ) effective 
chiral Lagrangian in the notation of Gasser and Leutwyler The constants 
L$ and Lio are known from the comparison between x?T and low energy hadron 
phenomenology. At the scale of the p meson mass they are L 8 = (0.9 ±0.3) x 10~ 3 
and Lio = (—5.5 ±0.7) x 10~ 3 . The high energy constants H\ and H 2 correspond 
to couplings which involve external source fields only and therefore can only be 
extracted from experiment given a prescription. 

3.3 The method and Ward identities 

The method used here is identical to the one used in ||. The full two-point 
functions are the sum of diagrams like those in figure [2]a. The one-loop two-point 
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(a) (b) 



Figure 2: The graphs contributing to the two point-functions in the large N c 
limit, a) The class of all strings of constituent quark loops. The four-fermion 
vertices are either £njl or ^njl i n ec L- ( P-2|) . The crosses at both ends are the 
insertion of the external sources, b) The one-loop case. 



functions are those obtained by the graph in figure Using a recursion formula 
that relates the n-loop graph to a product of the one- loop and the (n-l)-loop 
graph and the relevant combination of kinematic factors and Gy and G$ the 
whole class of graphs can be easily summed. Some care must be taken in the case 
where different two-point functions can mix so a matrix inversion is necessary 
(see ref. ||). 

The two-point functions defined above satisfy the following Ward identities. 
(We suppress the argument Q 2 for brevity.) 

-Q 2 n^ = (>/>,- >/>,)II^V (3.26) 

-g 2 n^. = ( mi -m j )U Sij + {q i q i )-{q j q j ), (3.27) 

-Q 2 n$\. = K + m,)n^., (3.28) 

-Q 2 IT^. = K + m J )n PiJ + (g 4 g i ) + (g J g i ). (3.29) 

These are derived in the appendix [B]. From these the flavour symmetry of the 
mixed two-point functions can be derived from the vector ones. 

The one-loop expressions, which we shall denote by II and use further the 
same conventions as given for the full ones above are given in appendix y. They 
satisfy the same identities but with the current quark masses rrii replaced by the 
constituent ones, Mj. In addition to these, there are two more relations that 
follow in general if the one-loop part can be described by a heat-kernel expansion 
in terms of the quantities E and R^ u of appendix [B|. These identities are (with 
the flavour subscript ij and argument suppressed) 

ny+nP = ng ) +ng ) , (3.30) 



8 



n 5 + g 2 4 0) = n P + Q 2 uf. (3.31) 



3.4 The transverse vector sector 

We introduce here for convenience an extra symbol gy 



9 v = j^fQ* . (3.32) 
The full resummed transverse vector two-point function is then 

rrCD 

n (l) _ LL Vij ( o oo) 

1 + gvU-vij 

This can be simply written in a form resembling the one in the complete VMD 
limit with couplings fg, fy and My depending on Q 2 and flavour and defined by 

rtW(r>2\ - 2 fs(Q 2 )jj 1 2 fy(Q 2 )ij M v(Q 2 )ij ( ~ n A \ 

" ~~Q^ + M 2 (Q% + Q 2 ' (3J4) 



2 fs(Q'% = ~ -(o . (3-35) 



• 2 , / y » , Q II y (Q j , , 

iV c A^ 1 
8ir 2 G v 1 . 

p2/^2\ _ tt(° +1 )^2 



TT ,! 7^2 

2/, 2 -W 2 )« = HT" (<?")« • (3-37) 



Where we have used the fact that (see appendices g and D) n^ +1) = U { y ] + ITy 
has no pole at Q 2 = 0. There is a correction here (in Ily ) due to the mixing 
with the scalar sector, which is allowed by the presence of explicit breaking of the 
vector symmetry (see the scalar mixed sector subsection |3.7] ). For the diagonal 

case, this is defined as = rrij or Mj = Mj, Hy vanishes and the formulas 
above simplify very much. 

The pole mass of the vector corresponds to the pole in this two point function 
or to the solution of Re (Q 2 + M 2 (Q 2 )ij) = 0. Alternatively, one can define the 
VMD values for the vector parameters (fy and My) as the best parameters of 
a linear fit of the inverse of U^\Q 2 ) i:i - 2fl(Q 2 ) ij /Q 2 . These definitions have 
the advantage that they are also valid for the Euclidean region (Q 2 > 0) where 
the vector cannot decay into two constituent quarks. See sections on numerical 
applications |3.9| and Vector- Meson- Dominance [5] for further comments. 
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3.5 The transverse axial-vector sector 

The transverse axial-vector two-point function derivation is also identical to the 
one in ref. ||. 

n i\ = -%) • ( 3 - 38 ) 

1 + gv^-Aij 



Using the identity (|3.30| ) it can be seen that this has a pole at Q =0 because 



has it. As can be seen from the explicit expression and is proved in general 
in appendix [B[ the combination Hy^ + Hy is regular at Q 2 going to zero. This 
again allows us to separate the pole at Q 2 = in a simple fashion. 

UAiQh ~ ~Q 2 ~ + M 2 A (Q 2 ) tJ + Q 2 ' (3 - 39) 
fUQ 2 ) = 9a{Q%JUQ 2 )i (3-40) 
VUQ 2 ) = -Q 2 U { a\Q%, (3.41) 



(9A(Q 2 h) = l-gvnf(Q% , (3.42) 



N A 2 

2f 2 AQ%M 2 A {Q% = ^^9A(Q 2 h, (3.43) 

o7T Kjy 

Va(Q% = g 2 A (Q 2 hn ( v +1 \Q 2 h ■ (3-44) 

There is a correction here (in H A ) due to the mixing with the pseudo-scalar 
sector due to the presence of both spontaneous and explicit breaking of the axial- 
vector symmetry (see the pseudo-scalar mixed sector subsection). For further 
discussion of these expressions and the ones in the previous section we refer to 



the subsection 3.8 on Weinberg Sum Rules. 



3.6 The pseudo-scalar mixed sector 

The same method as used in [H still applies with the results for the summed 
functions given in terms of the function Ap(Q 2 ) and the one loop two-point 
functions (with flavour subscripts ij suppressed), 



nf(Q 2 ) 

iW) 
rip(Q 2 ) 

A P (Q 2 ) 



A P (Q 2 



1 - g s U P (Q 2 ))U^'(Q 2 ) + g s (Up (Q 2 )) 



A P (Q 2 
1 

A P (Q 2 



1 <(Q 2 ) 



1 - g v uf (Q 2 ))U P (Q 2 ) + g v {jC (Q 2 )) 



(3.45) 
(3.46) 
(3.47) 



1 - g v uf (Q 2 ) ) ( 1 - g s Tlp(Q 2 ) ) - 9s9v ( lip (Q 2 ) ) " . (3.48 
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Using the identities for the one-loop case it can be shown that the resummed 
ones satisfy the Ward identities of appendix || with the current quark masses. 
To show this it is also necessary to use the Schwinger-Dyson equation for the 
constituent quark masses in eq. ( |2.3|) . 

In order to rewrite this in terms of a nicer notation we first express Ap(Q 2 )ij 
in a different form using the identities for the one-loop two-point functions. 

Ap(Q%- = "m'IV-' H {Q2) + Q2 ) (3 ' 49) 
withmJ(Q 2 ) = ( ";+m;L • (3-50) 



g S 9A(QWp(Q 



2\ 



I'.) 



Inserting the definition of f1j(Q 2 ) and 1/gs = — (QiQi) I '{Mi ~ m i) we recover 
the Gell-Mann-Oakes-Renner (GMOR) relation for the pion mass |T7| when eq. 
( |3.50| ) is expanded in powers of m,. For further discussion on corrections to the 



GMOR relation in this model we refer to the section on numerical applications [3.9| . 
Formula ( |3.50| ) gives the expression for the pole due to the lightest pseudoscalar 
mesons in the presence of explicit chiral symmetry breaking. 

This then allows us to rewrite the full two-point functions in a very simple 
fashion: 

ng>(q% = 2fi, { Q> ) ( mUQ l ) + Q2 -± s ), (3.5i) 

M 2 Mi+Mj 1 

AQh = 9s mim + Q 2 ' (3 ' 52) 

\2 



1 (Mi + MjY 1 1 

Ur(Q)ij ~ ~7s + Vim 7 s mUQ 2 ) + Q 2 - ( ' J, " J ' 1) 



Here we want to point out that the two-point functions lip 7 and lip suffer 
from the same ambiguity (via its dependence on g$) as the quark-antiquark one 
point-function (see discussion at the end of section ^) when compared with the 
xPT results. 

3.7 The scalar mixed sector 

Here we have to extend the analysis of to include possible mixing effects. This 
can be done in the same way as in the previous subsection with the result (with 
flavour subscripts ij suppressed), 



4 0) (Q 2 ) 



A 5 (Q 2 



1 - g s U s (Q 2 ))U^(Q 2 ) + g 8 (U$ (Q 2 )) 2 



(3.54) 



nf(Q 2 ) = x^W) n s {Q2h (3 - 55) 
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n s (Q 2 ) = 1 



(l -^< ) (g 2 ))n 5 (Q 2 ) + gv (nf (q 2 )) 



, (3.56) 



A 5 (g 2 

A 5 (g 2 ) = (i-^ni / 0) (g 2 ))(i-^n 5 (g 2 ))-^(nf(g 2 )) 2 .(3.57) 

To rewrite this in a simple fashion we would again like to expand As in a simple 
pole like fashion. Using the identities for the one-loop two-point functions this 
can almost be done, we obtain 



a*(q% = ys ^l y ((M, + Mj) 2 + UQ 2 W^ 2 ) + g 2 ) 



^np(Q 2 ^ 

— ^JWj f IV. 

Vdi — (71; 



+ n-(g 2 ) l ,^ 5 -^ I |— . (3.58) 

It can be seen that in the diagonal case a simple expression for the scalar meson 
pole can be found, 

Mj(-M|) _ = (Mi + Mjf + g A (-M 2 s ) lim U-M 2 s ) . (3.59) 
The expression for the scalar two-point function ITs(g 2 ) is in this case 

So in the diagonal case a simple relation between the scalar mass, the constituent 
masses and the pseudoscalar mass remains valid to all orders in the masses. In 
this case nj? = Ilf = 0. 

For the off-diagonal case, i.e. rrii ^ rrij, the corresponding expressions 
II^ and lis can be obtained from eqs. (|3.54j) - (|3.57|) and the explicit II functions 
in appendix There is a small shift in the pole compared to eq. ( |3.59| ) for the 

case rrii ^ rrij. From appendix |C], in eq. (|U.4j) , it can be seen that ITy itself has 

a zero close to a value of Q 2 = M| of eq. ( |3.59| ). In addition Uy^ is suppressed 
by (Mi — Mj) 2 /Q 2 . Therefore the value of the pole in the off-diagonal case is not 
too far from that in eq. ( |3.59|) . 



Here we want to point out that (as in the mixed pseudoscalar sector) the 
two-point functions Hy\ Tl<f and Us suffer from the same ambiguity (via its 
dependence on g$) as the quark-antiquark one point-function (see discussion at 
the end of section |2|) when compared with the %PT results. 

3.8 Weinberg Sum Rules 

The Weinberg Sum Rules are general restrictions on the short-distance behaviour 



of various two-point functions [IS|. They were first discussed within QCD in 
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ref. |T!J. A low-energy model of QCD should have a behaviour at intermediate 
energies that matches on reasonably well with the QCD behaviour. The general 
behaviour should be (JIlr = Ily — Ha-) 

lim (Q 2 nS°+ 1} (g 2 )) = First WSR, (3.61) 
lim (Q 4 nSj(Q 2 )) = Second WSR, (3.62) 
lim (Q 4 n£> (Q 2 )) = Third WSR . (3.63) 

Q 2 — *oo ^ ' 

Let us review first the QCD behaviour of these Sum Rules. In the large N c limit 
the three WSRs are theorems of QCD in the chiral limit (i.e., M. — > 0). The 
first WSR is still fulfilled in the large iV c limit with non-vanishing current quark 
masses. However the second and the third ones are violated as follows p0| , 

Jim^ (Q 4 II« (Q 2 )) = - hmjQ^UQ 2 )) 

= 2(m i (q j q j ) +m j (q i q i )) . (3.64) 

As shown in || the class of ENJL-like models does satisfy the three WSRs 
in the chiral limit. We shall now check how well this does in the case of explicit 
breaking of chiral symmetry. 

The high-energy behaviour of the two-point functions U-va needed for the 



three WSRs can be easily obtained from the expressions in sections |3]J, [3.5| , pi) 



and |3.7| . The first and second WSRs are satisfied in these ENJL-like models 
even with non-vanishing and all different current quark-masses. The high energy 
behaviour (Q 4 ) of these models is thus too strongly suppressed for ITjj^ (Q 2 ) to 
reproduce the QCD behaviour in the second WSR. The third one is violated as 
in QCD and one has 

lim (Q 4 n^(Q 2 )) = — (rriiMj + rrijM,) . (3.65) 

Q 2 ^oo v 9s 

Let us now see what relations between low-energy hadronic couplings do these 
Sum Rules imply for this ENJL cut-off model. In the equal mass sector, mi = 
rrij ^ 0, one has 

/ 2 M 2 = f\M\ + fl , (3.66) 
f v M 4 = f A M\ . (3.67) 

Remember that in QCD one has in this case 

f Y Ml = f\M\ + fl , (3.68) 
f v M v = f\M\ + m\ fl . (3.69) 

In the off-diagonal case, rrii ^ rrij, the situation becomes a lot more complicated. 
However, since the off-diagonal part is suppressed by (Mj — Mj) 2 /Q 2 one does 
not expect qualitatively different results. 
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Figure 3: The inverse of the transverse vector two-point function for equal quark 
masses in the chiral limit, i.e. M. — > 0; for the p meson, i.e. m\ = m<i = 3.2 MeV 
and for the <fr meson, i.e. mi = m2 = 83 MeV. The units of q 2 are GeV 2 



3.9 Some numerical results 

As can be seen from the explicit formulas the change with respect to ref. || 
is in most cases a (small) shift in the two-point function mass pole positions. 
Therefore we do not plot too many of the two-point functions. As numerical 
input we use for Gs, Gy and A x the values from fit 1 in ref. [Q]. These are 
A x = 1.160 GeV and G s = 1.216. The value of g A {Q 2 = 0) there was 0.61. This 
is Gy = 1.263. For the current quark masses we use the value of the quark mass 
for m = m u = that reproduces the physical neutral pion and kaon masses. 
With the other parameters as fixed above this ism = 3.2 MeV and m s /m = 26. 

As an example we have plotted the inverse of the transverse vector two-point 
function in eq. ( |3.34|) in figure ||| for the values of Gs and A x corresponding above 
mentioned. The full curve is the result in the chiral limit (M. — > 0) and the 



dashed is the result with m,- 



m,- 



m the value above. The reason we have 



plotted the inverse will become clear in section [5| We also show the inverse for 
rrii = rrij = m s the value above in the short-dashed curve. To show the result for 
unequal quark masses we have plotted in figure |] the transverse vector two-point 
function itself for the chiral limit case and for the us case with m s and m above. 
Notice that the two-point function now has a kinematical pole at q 2 = 0. 

We have also plotted in figure [| for the parameters quoted above the depen- 
dence of the pion mass on Q 2 . Since fijfn 2 ^ is a constant, see eq. ( |3.50|) this is 
also the Q 2 dependence of the inverse of the decay constant squared. 



14 




- 



-0.05 - ', : 

-0.1 - i ; 

-0.15 - | 

-0.2 I ' 1 1 1 h 1 1 

-0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 

q A 2 

Figure 4: The transverse vector-two-point function for the chiral limit and for 
unequal quark masses, m\ = m and m 2 = m s . Note the kinematical pole at 
q 2 = 0. The units of q 2 are GeV 2 . 
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Figure 5: The running pseudoscalar mass squared, m 2 -{—q 2 ), as a function of q 2 
for 7T7.j = rrij — 3.2 MeV. 
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Let us make some comments on the corrections we find to the GMOR relation 
( |3.50|) in this model. The corrections to the GMOR relation |l7j can be calculated 



here in an analogous expansion to the one in x?T. Then the GMOR relation can 
be written as follows [16| (for the diagonal flavour case, i.e. i = j) 



(3.70) 



2m,faft) = -mU-q 2 )fl(-q 2 ) (l - 4^M (2 L 8 - H 2 ) + 0(p G 



Here we have included all the chiral corrections to the quark condensate, to the 
pion mass and to the pion decay constant in their respective values. Then the 
remaining is a correction to the GMOR relation. We have also calculated this 
correction in this model and it turns out to be 

4^ (2L 8 - H 2 ) + 0(p«) = ^. (3.71) 

Notice that the r.h.s. contains all the orders in the x?T expansion in the large 
N c limit. Numerically, this correction is around 1 % for pions and 20 % for kaons 
and approximately agrees with the one found in QCD Sum Rules [2~I|. For the 
combination of 0(p 4 ) couplings 2L§ — H 2 in this model we get 

2Ls ~ H2 ~ i&? 2 r(-i, e )' (3 ' 72) 

with e = M 2 /A^ and M the constituent quark mass in the chiral limit. (For 
definitions of the incomplete Gamma functions V(n, e) see appendix 0.) The 
expression in ( |3.72| ) is the one consistent with the use of Ward identities to sum 
the infinite string of constituent quark bubbles. Numerically we get 2L 8 — H 2 — 
1.3 • 10~ 3 for the input parameters above. This differs from the one found at the 
one- loop level, in this same model in ref. (see footnote 1 ), numerically they 
find 2L 8 -H 2 = 0.2 • 10" 3 . 



4 Some three-point functions 

4.1 VPP with the use of the Ward identities 

In this subsection we calculate the Vector Pseudoscalar Pseudoscalar (VPP) 
three-point function to all orders in using the same type of methods as 

those used for the two-point functions. The three-point function we calculate is 
the following 

n^ PP (pi,p 2 ) = i 2 J d 4 x J d 4 ye i(pvx+P2 - y) (0\T (y^>(0)P kl (x)P mn (y)) |0) . (4.1) 
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Where k, /, m and n are flavour indices. In the limit of large N c the flavour 
structure is limited because of Zweig's rule (this flavour structure is general for 
any three-point function of three quark currents), 

n^ PP (Pl,P2) = ^(pl,P2)ikmSilS kn 5 mj + U~(p 1 ,p 2 )ik m SinSkjS m l . (4.2) 

Bose symmetry requires that 

n ^(Pl,P2)ifcm = n^(p 2 ,Pl)imfc • (4.3) 

The three-point function li^ pp (pi,P2) can then be simply calculated by only 
taking one particular flavour combination. Finally we can use Lorentz-invariance 
to rewrite 

U+{p u p 2 )ikm = Pl^fkmiP^Pl^ 2 ) + P2^fkm(pl,Pl,Q 2 ) , (4-4) 

where we have defined q = p\ + p 2 - 

We shall limit ourselves to the vector diagonal case, i.e. m; = m r In the 
vector off-diagonal case there will also be non-trivial mixings with the scalar- 
pseudoscalar-pseudoscalar three-point function. Here a relatively simple Ward 
identity for this three-point function can be derived from d^Vff = and the 
equal-time commutation relations. It is 

q'"U+ip 1 ,p2)ik m = -Up(-pl) ki + Tip{-pl)mk ■ (4.5) 

So the Ward identity relates the three-point function to a combination of two- 
point functions. This determines one of the two functions IT 4 , II 5 in terms of the 
other. The Ward identity gives, for instance, the following constraint (for p\ = p\ 
and i = m) 

nf fci (pV,g 2 ) = -Uf kl (p 2 ,p 2 7 q 2 ) . (4.6) 

The type of graphs that need to be summed are depicted in figure |6|. Each 
of the three tails here is the diagram in figure with the same explanation as 
there. We have there depicted one particular flavour combination. This is the 
one that corresponds to the function 11+ given above. The i, k, m written above 
the lines are the flavours of each line. 

All graphs are formed by having the tails summed over 0, 1, 2, • ■ -, oo loops 
connected by four-fermion couplings. These then couple to the one-loop three- 
point function (or vertex) 11^ , with various possibilities for the insertion in the 
three-point vertex. These possibilities for the 7-matrices are written in figure || 
inside the main loop. 

In this figure the left-hand side depicts the insertion of the current V^(0) and 
Tail I is the connection to this current. On the end connecting to the one-loop 
three-point function it is only nonzero for another vector insertion since in the 
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Figure 6: The graphs that need to be summed in the large N c limit for the 
Vector-Pseudoscalar-Pseudoscalar three-point function. See text for explanation. 



diagonal case we consider, the mixed vector-scalar two-point function vanishes. 
It expression is given by 



z 8fGv 



(4.7) 



Here the first term comes from where the external current directly connects to the 
one-loop three-point function and the second term is with the two-point function 
in between. The sum of both is 



~Q )mi 



(4.8) 



A similar discussion can be done for Tail II and Tail III. First we have the 
insertion of the current P kl (x) at the external end. On the end connecting to 
the one-loop three-point function we can have 275 or an axial-vector insertion 
since the mixed axial- vector-pseudoscalar two-point function is nonzero. The 275 
insertion tail is : 



47T 2 Gg 2 

1 + „ n P (-p 1 ) 



N c Al 



ki 



(M k + Mi 



29sfU-pi)(rnU-pi)-pi) ■ 



(4.9) 
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For the connection with the axial-vector insertion it is instead 

87l2Gv iy a n M ( v 2 ), 

N A2 Wl^P { Pl)ki 

^ X 



zp? (M fc + M,) 



2/£M£< 7 < ? (m 2 fei (-p 2 )-p 2 ) 



The combination /y M 2 is here flavour and p\ independent, it is the combination 
in eq. fl3.36| ) with Hy\Q 2 )ij = since we are in the diagonal flavour case. The 



way both these types of insertions can appear due to the tail are how within this 
formulation the mixing of pseudoscalar and axial- vector degrees comes about. 
These will be described by factors of g\ (see below). Tail III is identical to Tail 
II with the substitutions p\ —>■ p 2 and i, k — > k, m. 
The full expression for 11+ is 



8tt 2 GV y 2 
NcA l > 

X \K(Pl,P2) ( 1 + %^H P (-Pl) k ] U + ^n P (-pl) mk 



Y=vpa, f An 2 G s „ , 2 \ \ [8tc 2 G v . p M 2 

{pi,P2) \ i + -j^^Upi-p^ki I l ^ A2 zpgn P (-p 2 ) mfe 

X / X 



+nrr(p 1 ,ft) (^«n?(-ri)«) (i + igjf M-riw) 

(4.11) 

y pj± VAP VAA 

Where the one-loop three-point functions II , II and II are the one 
fermion-loop result for 

n^(pi,p 2 ) = i 2 J d 4 x J dSe^ + ^y\0\T(v^(0)P k \x)A^(y)) |0> , 

(4.12) 

n^ P (pi,p 2 ) = ? J d 4 x J d 4 ye i ^ x+p ^{0\T(v; m {0)A k j(x)P mk {y)) |0> , 

(4.13) 

n^/(pi,p 2 ) = * 2 / d 4 x / dV (pi - x+P2 ^(o|r(^r(o)^(x)Ar(i/)) |o>. 

(4.14) 



To obtain the full expression in eq. ( 4.11|) it now remains to calculate these 



VPP, VAP, VPA and VAA one-loop three-point functions (or vertices). The 
axial-vector ones always come multiplied with the relevant momentum. So we 
always have the scalar products p\ ■ A kl (x) and p 2 • A mk (y). That means that 
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using the Ward identities we can relate the VAA, VAP, VPA to the VPP one 
plus possibly two-point function terms resulting from equal time commutators. 
These Ward identities are (remember we assume Mi = Mj here). 

irfn™(pi,P2) = -(M k + M i )Tf^ A (p 1 ,p 2 ) 

-n Ma (-?)mi + n^t-paW ; (4.15) 
iplulf{ Pl ,p 2 ) = -(M k + Mi)U+( Pl ,p 2 ) 

+iU Pfl (-p 2 ) mk . (4.16) 

The other needed ones can be derived from this using Bose-symmetry. Notice 
that there is no contribution here from the flavour chiral anomaly (see eq. ( |4.24| )). 

We can now use these identities to obtain the final result for the three-point 
function we want. The terms which after the use of the one- loop identities above 
are proportional to VPP can be combined into a simple form using 5uG° 2 )- The 
result is (we have Mi = Mj and j = m in this flavour configuration). 



n + ^i,p 2 ) 



(Mi + M k f \ (g^M*(-q 2 ) mt - q»q 



A9 2 sfl(-pl)fL(-P 2 2)J V MU-q 2 U-q 2 J 

9A(-p 2 l)kigA(-pl)mkTC(PuP2) 



{m 2 k i(-pl) -p\) (m 2 mk (-pj) -pi) 

, (1 ~ 9A{-P\)ki) (1 - 9A{-Pl)mk) u , ( v X jdl) ( 2- 

+ 7T} , 7TT2 {(P2- q)Piu- (pi • q)P2u}^ v (~9 , 

{Mi + M k ) 

g A {-pl)ki{l ~ 9A(-Pl)mk) W M, 2 ^ 

Pl^llp {-PiJki 



mi 



Mi + M, 



k 



, gAj-pPmk (1 - 9A(-Pl)ki) =M. 2x 

1 ^r~T7 P2v^ P {-p 2 > 



mk 



Mi + M k 

(4.17) 

This result satisfies the Ward identity ( |4.5| ) if the one-loop function II one satis- 
fies the same one with the one-loop functions. This provides a rather non-trivial 
check on the result ( |4.17| ). 

It now only remains to calculate the one- loop form factor H (p\,p 2 ). We 
give its expression in appendix |B[ At this point we can see in eq. (|4.17 ) how 



far regularization ambiguities affect the result. We first have to define the two- 
point functions. Here all ambiguities are restricted to two bare functions (see 
section |3| for details). This three-point function adds one more in general, the 
three-propagator function h(p 2 ,P 2 ,Q 2 ) (see explicit expression in appendix 0). 
Of course, this one-loop form factor n^(pi,p 2 ), satisfies all the identities eqs. 
(|4. ip to ( |4.6p as well. We refer to section |5.2| for the definition of the physical 
vector form factor after reducing this VPP three-point function. We shall also 
discuss there the VMD limit in this form factor and give some numerics. 
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The same three-point function can be calculated in Chiral Perturbation The- 
ory. The result is 



1,P2 



2B ofn 



2L, 



(p2- P iT[l + ^q 2 + 0(p' 



{™>li ~ Pl)( m mk ~ Vl) 



f, 



2 

mi 



(4.18) 



Pulling out the pion poles (see section |5.2j for technical details) and taking the 
low-energy limit and the value of Lg in this class of models our full result in eq. 
(|4.17|) reduces to this, providing one more non-trivial check. 



4.2 PVV with a discussion about its Ward identity 

In this subsection we calculate the Pseudoscalar Vector Vector (PVV) three-point 
function to all orders in %PT with the same method as the one used before. 

n^ v ( Pll p 2 )=i 2 [ d 4 x J dV (pi X+P2 y) (0|T(P^(0)^(x)V7 m (y)) |0> . (4.19) 

Where i, j, k, I, m and n are flavour indices. A similar discussion about the struc- 
ture due to Zweig's rule can be given as was done before. We do the analogous 
decomposition into 11+ and H~ v functions as was done for the three-point func- 
tion VPP (see previous section). We shall here restrict ourselves to the case 
where all current masses or constituent masses are equal. Our main aim in this 
subsection is to show how the flavour anomaly affects the use of the one-loop 
identities. 

The class of graphs needed here is shown in figure [7] for the II + flavour com- 
bination. Each of the three tails here is the diagram in figure |2]a with the same 
explanation as there. The vector-like tails, II and III, can be easily summed (see 
discussion for the summation of the vector tail in the previous section) to obtain 
the overall factors 

X 

and 

9* + ^^Ki-PD ■ (4-21) 

Tail I can again couple at the one-loop end to both an axial-vector and pseu- 
doscalar two-point function. These have the same form as equations ( [4.9| ) and 
(|4.10|) in the previous section with pi — > q. Summing up the three tails the total 



result is then 
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Figure 7: The graphs that need to be summed in the large N c limit for the 
Pseudoscalar- Vector- Vector three-point function. See text for explanation. 



' r ll r ""( £)) (g^ + ^f^n v ^(-p 2 2 ) 



X 



An 2 G 



N c Al 

,P2 




tAVV 



(4.22) 



with II the one-loop result for the following three-point function 



n 



AW 



(Pl,P2) 



d A x \ d A ye 



i(pi-x+p2-y) , 



(0\T(A*™(0)V*\x)V u mk (y))\0). 



mk I 



(4.23) 



The main new part here is that at the one-loop level we now have to in- 
clude the anomalous part of the Ward identities. There has been in fact quite 
some confusion whether this can be done consistently. We have shown how this 
subtraction needs to be done in the case of ENJL-like models in ref. |J. The 
anomaly itself in this class of models is not well defined but a consistent subtrac- 
tion procedure to obtain the QCD flavour anomaly can be easily formulated, see 
ref. H . We could in principle use the prescription of ref. directly to obtain the 
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one-loop three-point function. Here we want to apply it to the PVV three-point 
function to all orders in external momenta and quark masses. The prescription is 
essentially to use the anomalous QCD Ward identities for the axial current con- 
sistently. We shall use the scheme where vector currents are conserved |23| . The 
subtractions these Ward identities impose in order to reproduce the correct QCD 
flavour anomaly in the ENJL-like models effective action we are working with, 
lead to the use of the following consistent one-loop anomalous Ward identity 

+ ^aUx)d a ai(x) + \i {v%(x), a% m ai k {x)} 
+ ^i{a a v^a p ) u {x) + ^{aWa^uix) 



with 



= 9V-9V-i[/,/] and 
dV = & i a v -zhAal . 



(4.24) 



Where and are the external vector and axial-vector sources defined in eq. 
( |2.1|) . As shown in ref. |§, when using this Ward identity, the fact that the 
anomalous part in eq. ( f4.24|) only contains external fields amounts to keeping 
the usual Wess-Zumino term |^4| (the only one of 0(p 4 ) in the chiral counting) 
for couplings of pseudoscalar type via G$ to external fields but when there are 
couplings to spin-1 fields via the Gy term, only the local chiral invariant part of 
the full term remains. We have checked that the form of the action given in M 
yields the same result as the one given below. 



So when we use the one-loop anomalous Ward identity in eq. ([4.24|) to reduce 
the right-hand side of Tail I to a part with only pseudoscalar couplings to the 
one-loop vertex, we obtain a local chiral invariant result plus an extra part where 
Tail I couples directly to the external vector sources {x)v™ n (y). This extra 
part is of order p A and is the subtraction the anomalous Ward identity imposes 
to obtain the correct QCD flavour anomaly. 

The full result in terms of the one-loop H three-point function is given by 




■ p x J [M v {-p 2 )u -p 2 
x <jl + ^_^n P (-g)-^f2M,n^(-g)^ 

+ XM,P2) ^^2M l U A P I (-q). (4.25) 
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Where the one-constituent quark loop function U is given by 



with F(plplq 2 ) = 1 + I 3 (plplq 2 )-I 3 (0,0,0) (4.26) 

where the form factor Iz{p\,p\, q 2 ) is the one given in appendix |B] and which 
appeared before in the study of the VPP three-point function in section (|4.1|) . 
This form factor coincides with the one found in the context of constituent quark- 
models (see for instance |p5[j ) when the cut-off A x is sent to oo. Here, this is 
a physical scale of the order of the spontaneous symmetry breaking scale and 
therefore we have to keep it finite. The anomalous Ward identities in eq. ( |4.24| ) 



is telling us that terms which are of chiral counting different to 0(p A ) have to be 
local chiral invariant 15) but they do not fix the regularization for those terms. 
We therefore use here consistently the same regularization for them as in the non- 
anomalous sector. At 0(p A ) the chiral anomaly also uniquely fixes the one-loop 
constituent chiral quark anomalous form factor to be the one in eq. (f4.26|) when 

Pl=pt = q 2 = 0m- 



Here we have used the anomalous Ward identity in eq. ([4 . 2 2|) . A naive use of 
the two-point functions and Ward identities would have led only to the first term 
in the sum in eq. ( |4.25| ) . The second term is the result of enforcing the validity 
of the QCD flavour anomaly. Substituting the results on the two-point functions 
in section |3| we can write down the following explicit expression 

H+ (r> n)- Nc r rPrf ^ m 



9A(~q 2 ) 



jsfli-q 2 ) {rnU-q 2 ) 



i F(r ? r? n 2) My(—pl)uM 2 (—pl)i 



(M 2 (-plh-pfj (m 2 (-pI)u-pI) 

(4.27) 

We refer to section f>.3\ for the definition of the physical anomalous 7r°7*7* form 



factor after reducing this PVV three-point function. We shall also discuss there 
on the VMD limit in this process and give some numerics. 

5 Meson- Dominance 
5.1 Two-point functions 

Here we shall discuss the vector case, the axial-vector case is similar. The trans- 



rrij = m 



.i 



verse vector two-point function in eq. ( |3.34j ) reduces in the diagonal case, 
(the off-diagonal case can be done analogously) to the following simple expression 
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with 2f 2 (-q 2 )M 2 (-q 2 ) = (5.2) 

o7T Cry 

and 2/ y (-g 2 ) = Tl { y\-q 2 ) . (5.3) 

In the complete VMD limit this two-point function has the same form but with 
fv and My constants. Let us see how complete VMD works in this model. For 
that, we shall study the inverse of Uy\—q 2 ), which in the complete VMD limit is 
a straight line. This function was plotted in section |3.9| in figure [3]. There we can 
sec that IT^ — q 2 ) in this model is very near of reproducing the complete VMD 
linear form. Moreover, we can perform a linear fit to the inverse of H ( y\-q 2 ) to 
obtain the best VMD values for the fy and My parameters. These parameters 
are in this way meaningfully defined in the Euclidean region — q 2 > where the 
model is far from the two constituent quark threshold. Doing this type of fit for 



the values of the input parameters A x , Gy, Gs discussed in section leads to 
My ~ 0.644 GeV for the vector mass in the chiral limit (remember that we are 
always in the large N c limit) and fy ~ 0.17 for the decay constant. For current 
quark masses values discussed also in section |3T9| , we obtain for the p meson 
flavour configuration M p ~ 0.655 GeV and f p ~ 0.17 and for the <j) meson one 
Mtj, ~ 0.790 GeV and f^ ~ 0.14. We see thus that the p mass is very close in the 
large N c limit, to the one in the chiral limit, My. Notice that these values for 
My are far away from those quoted in ref. The underlying reason is that in 
ref. Q fy and My were determined directly from the Lagrangian at 0(p 2 ) in the 
ENJL expansion, identifying them with their values at q 2 = 0. What we find here 
is that even though the two-point function in eq. ( |5.1| ) has the correct q 2 — > 
limit behaviour it does have, with the choice of vector fields to represent vector 
particles in ref. M , substantial contributions from higher order terms (mainly of 
0(p 4 ) in the ENJL expansion). A physical vector field that would include these 
contributions can in principle be defined as is shown by the fact that the inverse 
of Hy (— q 2 ) is a rather straight line. What has happened is that 

f 2fyMy) 

M v (0) -q 2 (l + X + 0(q 2 /A 2 x )) • 

The vector meson mass derived in was My(0) while the slope of the physical 
two-point function (for |g 2 |/A 2 << 1 that is where this ENJL cut-off model makes 
sense) corresponds to rather My ~ M v (0)/\/l + A. We find from the calculation 
that indeed A is of order 1 (A ~ 0.7), explaining the difference in the slope from 
the 0(p 2 ) ENJL calculation in ref. 0] of the two-point function to the 0(p A ) one. 
We can also see from eqs. ( |3.39| ), ( |3.5ip -( |333"D and ( |3.60| ) that the forms 



of these two-point functions are very similar to the corresponding ones in the 
meson dominance limit but with couplings varying with q 2 . The identification of 
the corresponding physical values will involve analogous procedures to the one 
described above for the transverse vector two-point function one. 
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5.2 VPP three-point function 



In this subsection we discuss how the result for the three-point function I\X pp {pi ) p 2 ) 



obtained in section [4J] can be used to determine the physical pion electromagnetic 
form factor in this model. We shall discuss the VPP three-point function flavour 
structure corresponding to the three-point function n+(pi,p 2 ) i n ec L- (14.171) for 



m = rrii = rrij = and p 2 = p\ = p\ for definiteness. 

Since this n+(p 1 ,p 2 ) is a Green's function we first have to reduce the external 
legs to properly normalized pion fields. The vector leg acts here as an external 
source and is properly reduced without bringing in any factor. For this, we first 
look at the pseudoscalar two-point function in eq. ( |3.53| ) obtained using the same 
external fields and parametrize it around the pole as 



n, 



-p 



1 

9s 



1 + 0{m\lk\ 



(5.5) 



The reducing factor Z v is 



Z^ 
A 2 



'Vli-Oal A 2 

dm 2 A-p 2 



with 



1 + 



dp 2 



-s2 



(5.6) 



where I 3 (pl,p 2 , q 2 ) defined in appendix [D] and f 2 (—q 2 ) and f^—q 2 ) in eqs. ( [3.401 )- 
( |3.41 ). The quantity A is very close to one and exactly one in the chiral limit. 
Each pion leg brings a factor Z^ 2 after reducing the Green's function to the 
physical amplitude. Rewriting the pseudoscalar two-point function in the form 
in eq. ( |5.5| ) gives that m 2 is the solution of m 2 = m 2 j(—m 2 ). 

Reducing the VPP three-point function U+(pi,p 2 ) in eq. ( [4.17|) we find that 
it can be written as followso (we shall suppress the flavour indices which are 
always ii) 

Z-,r 



IT^(pi,p 2 ) 



Fvpp(p ,q ) (P2 -Pi) 



(5.7) 



(P 2 - ml) 2 

which defines the electromagnetic pion form factor (or in general the pseudoscalar 
vector form factor) Fypp(p 2 ,q 2 ) in this model. (The general pion form factor, 
i.e different quark masses and p\ ^ p\ can be obtained similarly from ( [4.17|) .) 
This form factor in the ENJL model is expected to be a good approximation 



2 To obtain the 7*7r + 7r three-point function from this 11+ is necessary to multiply it by the 
electric charge of the pion. 
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Figure 8: The inverse of the vector form factor of the pion of eq. ( |5.8|) . For the 
chiral limit and with all current quark masses equal to 3.2 MeV. Also plotted is 
the VMD approximation My(—q 2 )/(My(—q 2 ) — q 2 ) for the latter case. 



at intermediate and low-energy energies, within the validity of the ENJL model 
we are working with, i.e. for \q 2 \ « A 2 . The explicit expression for this form 
factor! 



is 



F VPP (ml,q 2 ) 



x 



1 M 2 {-q 2 ) 
2A 2 f 2 (-ml)M 2 (q 2 )-q< 



q \l-g A {-ml)Yf v {-?) + 



q 2 — Am 2 



(q 2 - 2ml)(fl(-q 2 ) - /'(-to*)) - 4m£l 8 ("4, ™l f)] } ■ 

(5.8 



Notice that this form factor has no pole at q 2 = 4m 2 . The value of A 2 in eq. 
( |5.6|) is precisely the one that ensures that F VP p(m 2 , 0) = 1 in the large N c limit 
as is required by the electromagnetic gauge invariance. This must be so since we 
have imposed the Ward identities to obtain this form factor. In figure [8] we have 
plotted the inverse of this form factor for the parameters quoted in section |3.9| in 
the chiral case (m = 0) and in the case corresponding to the physical pion mass 
(to = 3.2 MeV). As can be seen from the picture, it is a rather straight line so 



3 This form factor was also calculated in rcf. [||. With the appropriate changes of notation 
it agrees with the one found there. 
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the complete VMD result for this form factor, i.e., 



M 

b V pp (m n , q ) - 2 - 2 [b.y) 
p 9 

with constant vector mass M p works rather well. The slope of the linear fit of 
the inverse of the form factor in eq. ( |5.8| ) to this VMD form gives a vector mass 
which is M p ~ 0.77 GeV. This mass is very close to the physical value and rather 
different from the one found for the transverse vector two-point function in the 
VMD limit M p ~ 0.655 GeV in the large N c limit. This explains why using the 
physical p meson mass and the VMD dominance works so well but it also shows 
that this M p "mass" in eq. (|5.9|) has not, in principle, to be the same as the mass 
of the vector meson described by the transverse two-point vector function. 

The same three-point function VPP also contains implicitly the p — > tttt 
coupling constant gv- (See ref. Q for its definition. Notice that is different from 
the symbol used in section |3|.) Again, to obtain the physical p — > tttt amplitude 
we should first reduce the vector leg that now corresponds to the p particle, 
(remember that the pion legs have been already reduced). This will bring a 
factor which is similar to the factor 1 + A discussed in the previous subsection. 
We shall, as before, first determine the reducing vector factor from the vector 
two-point function in eq. (|5.1|) . The reducing factor Z p is 



VyM 2 y 



q 2 =M$; 



(5.10) 



In this equation the combination 2f v M v is the one given in eq. ( |5.2| ) and is 
independent of q 2 . The vector mass M p is again given by the solution to M 2 



One also can rewrite down the electromagnetic pion form factor showing ex- 
plicitly the coupling constant of the p meson to pions, gv, as follows 

q 2 M 2 

F VPP = 1 + fv9v ~Pi M i P _ 2 ■ ( 5 - n ) 

Then, in the complete VMD limit one has fvgv — fn/M 2 - In this ENJL model 
this relation is equivalent to gv — (1 — 9a) fv, i>e. one has complete VMD and 
the KSRF relation [[□]] 2gv = fv satisfied for g& = 1/2. 

One can see in the eq. ( |5.11| ), that reducing the p vector leg brings in a 
factor B 2 in the numerator and another factor B 2 in the denominator with the 
net result that fv{.~1 2 )9v{.~1 2 ) is n °t affected by reducing of the vector leg as 
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much as happens to fy(—q 2 ) in eq. 
eq. ( p. HQ we get the following 



( |5.1Q . Then, with the definition of fv9v in 



fv(-q 2 )9v(-q 2 ) = 
1 



2A 2 q 2 



q 



M 2 ( 



(2A 2 f 2 (-ml)/M 2 (-q 2 ) - f 2 (-q 2 )(l - g A {-ml)) 2 ) 



+(l + g A (-ml))f*( 

x((q 2 -2ml)Ql(- 



-mi 



+ 

7l( 



q 



Ami 



q 



(5-12) 

Notice that this fv(—q 2 )gv(—q 2 ) form factor has neither a pole at q 2 = nor at 
J - 4m 2 and when expanded in q 2 and with m\ = one gets /y(0)gy(0) = 2L 9 , 
where Lg is the one found in ENJL in ref. [Q. As discussed there, at q 2 = one 
has the KSRF |TTj] relation, i.e. /y(0) = 2g v (0) (which is valid for q 2 = m 2 = 0) 
analytically for g A = and very approximately satisfied for g A varying between 
and 1. The expression in eq. (|5.12j) is the off-shell equivalent to the KSRF 
relation in this model. For g A = the vector mass vanishes and the p meson 
couples as an SU(3)y gauge boson, in fact in this limit one recovers the results 
of the Hidden Gauge Symmetry model |27| for the non-anomalous sector. In 
particular, when g A = we have that the reducing factor B is 1 as corresponds 
to external gauge sources. In this limit (g A = 0), one still has the KSRF relation 
analytically satisfied off-shell, i.e. fv(—q 2 ) = 2gv(— q 2 ) for all q 2 . 
In the limit g A — ► 1 one has 



fv(-q 2 )9v(-q' 



A 2 q 2 



/;(-<)(! -a 2 ) + - 

q 



1 



4mj 



x 



((q 2 - 2ml)(f 2 (-q 2 ) - f 2 (-ml)) - Am%(ml m 2 n , q 2 )) 



(5.13) 



This is the constituent quark model result (in g A = 1 the vector mesons decouple 
from this model) and when expanded in q 2 with m 2 = it coincides with the 
corresponding result in ref. [J3J. The KSRF relation is not analytically fulfilled 
in this limit but on can see that analytically is very approximately satisfied. 
Then, we have that for g A varying between (the gauge vector limit) and 1 (the 
constituent quark limit) the KSRF relation goes from being analytically fulfilled 
to be very approximately fulfilled for any value of q 2 . 

Let us see how gv{~q 2 ) works numerically compared with fv{~q 2 ) for a def- 
inite value of g A . In figure ^| we plot fv(—q 2 )/2 and gv(—q 2 ) for the values of 
parameters discussed in section |3]9|. These values correspond to g A (0) = 0.61. 
The form factor gv(~q 2 ) is somewhat dependent on q 2 with (2.1 ~ 2.2) gv{~q 2 ) — 
fv{—q 2 ) in the Euclidean region. In this figure we also plot the case g A — > 1 where 
the same features can be seen. The form factor gv{—q 2 ) for any value of g A will 
be between the line fv/2 (i.e., the g A = limit) and the line for g A = 1, therefore 
the KSRF relation is approximately satisfied off-shell for any value of g A . 
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Figure 9: The generalized KSRF relation. We plot gv(—Q 2 ) for = 0.61 (solid 
line); qa — ► 1 (short-dashed line) and /y(— <? 2 )/2 (dashed line). The difference 
between the curves gives the violation of the KSRF relation. See text for further 
comments. 



5.3 PVV three-point function 

In this subsection we want to study the 7r°7*7* anomalous form factor. For that 
we shall reduce the PVV Green's function in eq. Q4.27D calculated in section 



47| to the physical amplitude following the same procedure that in the previous 
section (for details see there). Now, we have to reduce one pion leg, this will 
bring in a factor y/Z w and two external vector sources legs which are properly 
reduced without bringing any factor. Then the PVV three-point function in eq. 
can be rewritten as follows 



n+(pi,p 2 ) = — —rr^ie^ P PiP P 2 



x F PVV (q 2 ,pipj) (5.14) 

where Fpyy is the ir° — > 7*7* form factor in this model. Notice that the reducing 
factor A in eq. ( |5.6|) goes to one in the chiral limit preserving, in that way, the 
chiral anomaly condition Fpw(0, 0, 0) = 1. This form factor can be used as 
an accurate interpolating expression in low-energy hadronic processes valid for 



4 To obtain the 7r°7*7* three-point function from this is necessary to multiply it by 
a factor \[2 coming from the n° flavour structure and a factor e 2 /3 from the quarks electric 
charge. 
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Figure 10: The inverse of the 7r°7*7 form factor for one photon on-shell 
and one off-shell as a function of the photon mass squared, q 2 . No- 
tice the linearity in the Euclidean region. Plotted are the full result, 
M^(-q 2 )/(M^(-q 2 ) - g 2 )(VMD-like) and the ENJL model without vector and 
axial- vector mesons (g A — !)■ 



external momenta smaller than A 2 . We plot the inverse of this form factor for the 



case p\ = in figure |10|. Notice that there Fpyvijn^, 0,0) ^ 1 and the difference 
comes from the reducing factor A and is of chiral counting 0(p 6 ). We can expand 
this form factor for small p\,p\ and pion mass i as follows 

F PVV (m 2 n ,pl,p 2 2 ) = l +p (p 2 l+p 2 2 ) +p ' m l + 0(q 4 ), (5.15) 

this expansion defines the slopes p and p' which in this model are 

1 r(2,M 2 /A 2 )\ 
P = 9A{U)\— 277TT + 1 



M^(0) 12M 2 J 

:t(2,m 2 /a 2 ) t(i,m 2 /a 2 ) \ 

and p' = ^(0)1-4^- ^;^/^ ). (5.16) 

Where the second term in p' comes from the reducing factor A defined. The 
constituent quark mass M here is the one corresponding to the current quark 
mass value m = 3.2 MeV used in the numerical applications section |3.9| . Using 
M£(0) = 6M 2 5u(0)/(l - g A (0)) we can write down them as 

P = 77^72 (2 " (2 - T(2, M 2 /A 2 )) g A 



12M 1 



5 For the 7r° decay we are on the pole and hence q 2 = m 2 



31 



and 



j = !M (Y ( 2 a/ 2 /a 2 ) - r(1 - M2/A " 



12M 2 V x r(o,M7Aj) y 

(5.17) 

which interpolate between the constituent quark-model result #a(0) = 1 and the 
gauge vector meson result <?a(0) = 0. 

With the input parameters we have been using (see numerical application 
section ^.9| ) we get 

p = (0.86 + 0.67) = 1.53 GeV" 2 
and p = (0.67- 0.27) = 0.40 GeV" 2 . (5.18) 

Where for p the first number between brackets is the vector meson exchange con- 
tribution and the second is the constituent quark contribution (up to ^(0)). We 
see that both contributions are very similar giving some kind of complementarity 
between both approaches and explaining the relative success of both when used 
to describe this slope. For p' they are the constituent quark contribution and 
the one coming from the pion leg reducing factor 1/A. (Notice the cancellation 
there.) 

In the limits g& — ► 1 and — ► we find 





P 


= 1.10 GeV" 


-2 


for 


9 A ~ 


-> 1 




P 


= 2.20 GeV~ 


-2 


for 


9 a - 


-> 




i 

P 


= 0.66 GeV~ 


-2 


for 


9 A ~ 


-> 1 


and 


P' 


= 0.00 GeV~ 


-2 


for 


9 a ~ 


-> 



(5.19) 

We see that the difference between these two limits is big and that the actual 
result is some kind of interpolation. Experimentally p8| 

p = (1.8 ±0.14) GeV" 2 . (5.20) 

Taking into account that the 1/N C corrections from x?T loops are estimated p9| 
to be twice the experimental error we consider the result as good. 

Let us compare this full result in eq. (|5.17|) with the one obtained in ref. [[30 



in this same model assuming complete VMD in the chiral limit. There, the same 
prescription to include the QCD chiral anomaly that here was used at the 
one-loop level with the result 

" = mp^Aor- (5 - 21) 

Of course, this complete VMD result vanishes when — 1 where vector mesons 
decouple. The differences between eq. ( 5.1 7]) and eq. Q5.21 ) come from the 



resummation of all the orders in external momenta that are included in the full 
result in eq. fl5.17T) . One can see that the complete VMD result coincides with 
the full result for (?a(0) ~ 0.50. 
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6 Conclusions 



In this paper we have derived in a general class of ENJL-like models the two- 
point functions away for the chiral limit and for different masses in terms of the 
one-loop ones. This derivation used the Ward identities of the one-loop functions. 
The heat kernel expansion yields two more identities then can be derived from the 
current identities directly. These are then used to rewrite all two-point functions 
in terms of two basic ones. These can then be calculated in the bare ENJL-model 
as we did here or gluonic background corrections can be taken into account (see 
ref. H for a discussion). These extra identities allow us to discuss the Weinberg 
Sum Rules in this class of models. We find that the first WSR is satisfied and 
the third one is broken in the same way as required in QCD. The second WSR is 
satisfied in this model while QCD requires it to be broken. The very high energy 
behaviour is thus a little too suppressed. 

We also find (in the case of equal masses) the simple formula for the scalar 
mass that had been argued before when only divergent terms were kept in the 
heat kernel expansion (see second paper in ref. 0). All two-point functions 
can also be written in a form very like a form of meson dominance but with the 
couplings and masses depending on the momentum. In the Euclidean region, 
q 2 < 0, the vector two-point function can also be well described by a VMD form 
with constant couplings. The relation of these with those from the low-energy 
expansion was treated as well. 

We then proceeded to calculate two examples of three-point functions. Again 
the use of Ward identities simplified the calculation and pinpoints all the reg- 
ularization ambiguities into the one-loop function. We want to point out that 
both the anomalous sector and the non-anomalous sectors of these ENJL-like 
models are then treated on the same foot and VMD can discussed in both sectors 
with a unique prescription, namely the use of the relevant Ward identities ||, U . 
Then the regularization dependence uncertainties are consistently treated within 
the same prescription in both the anomalous sector and the non-anomalous sec- 
tor. Here we discussed Vector Meson Dominance and the KSRF relation for the 
VPP case. We then use the Ward identities, modified to reproduce the QCD 
flavour anomaly, to calculate the PVV function. Here we find that naive VMD 
expectations for this function cannot be realized. No simple generalization to q 2 
dependent couplings is possible. Formally our expression looks very much like the 
VMD expression with couplings and vector mass running with q 2 times Qa{— Q 2 ) 
plus a second term coming from the requirements of the anomaly. The numerical 
result for the slope is, however, in good agreement with the VMD value. But the 
5u( — Q 2 ) factor diminished the "real" VMD part to a little more than half the 
full value while the constituent quark loop adds the remainder. Here we reconcile 



both explanations for the slope, the VMD one and the quark loop one, [25, 2£ 
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A Proper time and incomplete Gamma func- 
tions 

In this appendix we give the regularization method we have used throughout 
this work and some related definitions. After performing the standard Feynman 
parametrization, one constituent fermion propagator is regulated consistently in 
this work using a proper time regulator as follows 



f°° dre-™ 2 ^V). (A.r 
Ji/M 



M\Q\x) A/a 

After performing the remaining Q 2 integration and the change of variables rA 2 
z one arrives to the following type of integrals 



X 



r(n-2,e) = / — z n - 2 e-\ (A.2) 

with e = M 2 (Q 2 , x)/A 2 x and n = 1, 2, • ■ -. These T(m, y) are the so-called incom- 
plete Gamma functions. 

In general, this regulator breaks the Ward identities. We have, however, 
always imposed all the Ward identities explicitly so our results have the correct 
symmetry covariance. 



B Derivation of the Ward identities 

In this appendix we generalize the proof in the appendix of ref. to the case 
with nonzero current quark masses. There a proof was given of all relevant 
identities in terms of the heat kernel expansion (for an excellent recent review 
and definitions see ref. |!T|) and some of them in terms of the Ward identities as 
well. Here those which can be derived directly from the Ward identities can also 
be derived from the heat kernel expansion but since they involve different masses 
they require a resummation of different terms. For these the direct derivation of 
the Ward identities is actually simpler. Only for the additional relations will we 
give the heat kernel derivation. 

At the one-loop level we use as Lagrangian the one in eq. (|2.2|) (with the same 
definitions as there) 

C ENJL = qDq (B.l) 
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where T> contains the couplings to the external fields l^r^s and p as well as 
the effects of the four-quark terms in £ N 'J L and £nj L on the quark currents at 
the one-loop level. In particular it contains the constituent quark masses, Mj. 
However, we shall keep the notation l^r^, s and p to denote the quark current 
sources in the presence of these four-quark NJL operators. The one-loop current 
identities derived from this Lagrangian are 

d^A^ = (M; ■ P ij . (B.2) 

When the whole series of constituent quark bubbles are summed these identities 
are satisfied changing constituent quark masses by current quark masses. In 
addition we use the equal time commutation relations for fermions 

WO*), <$(!/)} =iS a p6i j 6 3 {x-y) . (B.3) 

Here a and (3 are Dirac indices and x means the spatial components of x. Mul- 
tiplying the two-point functions with iq^ is equivalent to taking a derivative of 
the exponential under the integrals in eqs. ( |3.5| ) to (|3.8|). By partial integration 
we then get several terms, those due to the time ordering which leads to equal 
time commutators and those where the derivative hits one of the currents. The 
first type are evaluated using eq. (P-3|) and the second type are related to other 



two-point functions using eq. ( |B.2| ). This then leads to the expressions (|3.26| ) to 

(EE27D- 



The derivation of the other two identities is slightly more complicated. The 
effective action of the Lagrangian in eq. ( |B.1| ) can be obtained in Euclidean 
space as a heat kernel expansion (see ref. ||31|| ). The coefficients of this expansion 
are the so-called Seeley-DeWit coefficients, they are constructed out of the two 
quantities E and R^ u . These are defined as 



ptp = _v M V M + £ + M 



■2 



= [V„,V„] , 

= a M #-zK,#]-z[a M 75,#]- (B.4) 
If in eq. ( |B.1|) the Dirac operator T> contains couplings to gluons these should not 



be taken into account in eq. (|B.4|) . The relevant heat kernel expansion in that 
case will have different coefficients depending on vacuum expectation values of 
gluonic operators, but will still be constructed out of the quantities in eq. ( |B.4| ) 
(depending now also on the gluon field). The quantity M is the mass that is used 
in the heat kernel expansion. The operator T> is 

i^{d^ - iVfj, - ia^ 5 ) - M - s + ip^ 5 . (B.5) 
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Here M. = diag(m u ,md,m s ) is the current quark mass matrix and we allow for 
spontaneous chiral symmetry breaking solution (0|s(x)|0) 7^ 0. For the terms 
relevant to two-point functions we have 



and 



- ^ {^75, M + s- i]Tf 5 } + {M, s} - i [M, p] 75 + s 2 + p 2 

+ M 2 — ~M 2 , 

with 

v w = 5V-3V-i[/,/], 

d^# = ^#-iK,#]. (B.6) 

The main difference with ref. is the occurrence of the last line in the ex- 
pression for E in (|B.6|) . We shall call this last line E . In this equation, 



M = diag(M U} Md, M s ), the diagonal matrix of the constituent quark masses 
defined in eq. ( |2.3|) . Notice that the scalar field here has been shifted and we 
have now (0|s(x)|0) = (though we use the same notation for it). When Gs — > 
in eq. ( [2.3|) then M — > and M — > M.. Let us now systematically go through 
all possible types of terms in the expansion. We shall not discuss the mixed 
two-point functions here since we only want to prove eqs. (B.30|) -( |3~3TD . 



In the heat kernel expansion, those terms containing two factors R^ v only 
contribute to the transverse parts, Hy\ and in the same way. Their contributions 



hence obviously satisfy eqs. (|3.30| )- (|3.31|) . Similarly, one factor R^ u requires the 



presence of two covariant derivatives V M . By commuting derivatives (the extra 
terms only contribute to three and higher point functions) and partial integration 
these can be brought next to each other so they convert into a second factor R^ u . 
This brings us back to the previous case. Intervening E's can only contribute 
via Eq but these do not spoil the above argument. The first term in E, namely 
a^Rv, requires a 2nd a^ u R^ v because otherwise the trace over Dirac indices 
vanishes. These also behave like terms with two factors R^ v . Therefore, in the 
remainder we are only concerned with E without this first term. 

E can also directly contribute to the scalar and pseudoscalar two-point func- 
tion in the same way via s 2 + p 2 . Extra factors E become again Eq and extra 
derivatives also respect the relation (|3.31|) . The most complicated case is where 



both fields come from a different E. This contributes in the form EqEEq 1 cP 1 E. 
These contribute to all form factors in the form M^M^q 21 times the coefficients 
listed in Table [TJ. These coefficients obviously satisfy the relations ( |3.30| )-( p.31| ). 



The last type of terms is where the external fields come out of a derivative. We do 
not consider the mixed case here, so both the fields have to come out of a deriva- 



36 



Function 


Contribution 


n 5 
n P 

TT {0) 

TTd) 
11 A 

n (o) 

n (D 


-q 2 + {Mi + Mj) 2 
-q 2 + (Mi - Mj) 2 

(Mi + M 3 ) 2 /q 2 
-(M. + Mrf/q 2 

(Mi - Mj) 2 /q 2 
— (Mi - M 3 ) 2 /q 2 



Table 1: The contribution of terms of the type E m+n+2 to the two-point functions. 

tive due to the 7 M that is necessarily present in the E that would be a candidate 
for the external field. So there are those where the external fields are contained in 
two factors V M . If the indices of these are different, then there need to be at least 
two extra derivatives present that will produce a q^qu- This contributes equally 

to 11^ and 11^ . If the indices are equal, it will contribute proportional to 

and thus to the vector and axial- vector equally with ITy^ 1 ' ) = 0. This completes 
the proof of the identities (CT)-(CT). 

Now it remains to prove that these contributions will never produce a pole 
in n^° +1 ' ) at q 2 = 0. Terms that contain two factors R^ u contain two factors 
of momenta and hence do not. Terms with one factor R^ v can be brought in 
the form with two so do not produce a pole either. From Table [l] there is no 
contribution from that type of terms to Ily^ 1 . Then those with external fields 
from with different derivatives necessarily contain extra factors q^ so do 
not contribute to a possible pole at q 2 = and the last type of terms does not 
contribute to n^ 1 ' ) as shown above. This completes the proof. 



C Explicit expressions for the barred two-point 
functions 

Here we shall give the one-constituent-quark-loop expression for the two-point 
functions defined in eqs. (|3.5|) - (|3.10|) in the presence of current quark masses. 
These two-point functions are denoted in the text as the II ones. They fulfil the 
same Ward identities as the full-ones in eqs. (|3.26| )- (|3.29|) changing the current 



quark masses there by the constituent quark ones. In addition, they also satisfy 
the Ward identities in eqs. ( |3.30|) - (|3.31|) . Using these identities one can see that 
there are only two independent functions out of U.y , ILy , n^, 11^, n^f n^f, 
Hs and lip. We shall take lip 1 and 11^ + 11^ as these functions. The explicit 
expressions are 
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\ N r 1 

(Q% = ttt ^8 / dxx(l-x)T(0, Xij ), (C.l) 
/ 1d7t Jo 

K(Q%j = -^ £ 2 4/ 1 dx(M l x + M,(l-x))r(0,x u ), (C.2) 

107T J 



167T 2 7o 

where 



Mfx + M?(i-g) + g 2 g(i-g) rr ,,v 

X 



One can obtain all the others one-loop two-point functions in function of these 
gets 



two by using the Ward identities mentioned above. For instance, for the Hy^ one 



(o) (Mi - M£ < (Q%- 

U V W)lj M . + M . Q2 {Q 2 + (M . _ M . )2) 

x {(Mi + + ^(O 2 )^^) (l - (^*) (^|)) + Q 2 } • 

(C.4) 



D Explicit expression for the one-loop form 
factor n+(pi,p 2 ) 

Here we shall give the one-constituent-quark-loop expression for the three-point 
function n^(pi,p 2 ) defined in eq. (|4.3|) . We shall give it for Mi = Mk = M m . 
The explicit expression is (remember that we have j = m), 



n "(pi,P2) 



2M 



n 



+ 



■Mi 
L P {' 

Pl -P2 



-Pih 



PlP2 - (Pi ■ P2 

+ (Pi ■ P2) (n. P (-pl)ii - Rp(-q 



+ -j^h{pl,pl,q 2 )pl 



pl ( n p {-pDu - n p (-<? 2 )« 
Pl + (pi ■ P2) 



1 + (pi -P2)- 



P1P2 - (Pi • Pa)' 



(Pi 



-P2) • 



>P? 



:D.r 
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Where the two-point function n p (— p 2 ) was given in appendix and the function 



M 2 (x,y) = 

M 2 - p\{\ -x)- plx(l -y) + (pi(l -x)- p 2 x{\ - y)f . 

(D.3) 
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